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A  Mixture  of  Elastic  Continua 


A. 


* 

E.  Green 


by 

and  P.  M.  Naghdi+ 


1.  Introduction 

Mixtures  of  chemically  inert  elastic  continua  have  been  discussed 
previously  by  a  number  of  writers.  It  is  the  purpose  of  the  present 
note  to  provide  a  fairly  simple  discussion  of  the  constitutive  equations 
for  a  mixture  of  v  elastic  continua  with  a  single  temperature. 

Green  and  Naghdi  [l]  have  recently  made  a  small  change  in  the  use 
of  their  original  theory  of  interacting  continua  [2]  involving  a  single 
temperature.  On  the  basis  of  this  theory,  we  consider  here  a  mixture  of 
v  chemically  inert  elastic  continua  whose  constitutive  equations  are 
nonlinear  functions  of  temperature  and  suitable  kinematical  variables  of 
each  constituent  but  are  linear  functions  of  degree  one  in  temperature 
gradient  and  velocity  differences.  The  results  obtained  have  a 
simple  form. 
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2.  Notation  and  preliminaries 


We  consider  a  mixture  of  v  interacting  constituents,  each  of  which 

is  regarded  as  a  .continuum.  We  assume  that  each  point  within  the  mixture 

is  occupied  simultaneously  by  all  v  constituents  which  are  °in  motion 

relative  to  a  fixed  system  of  rectangular  Cartesian  axes .  The  position 

till 

of  a  typical  particle  of  the  a  constituent  at  time  t  is  denoted  by 
x®(t)  ,  where 

xi(T)  “  Xi(xi»  *2*  x3>  t)  (-•  <  T  *  t)  ,  (2.1) 

a 

and  is  a  reference  position  of  the  particles  of  this  constituent. 

All  Latin  indices  are  tensor  indices  and  take  the  values  1,2,3  and  the 
usual  summation  convention  is  employed.  'Ihe  Greek  letter  a  ,  as  in 
(2.1),  is  reserved  for  reference  to  the  constituent  of  the  mixture 
only.  Unless  indicated  otherwise  the  summation  convention  does  not  apply 
to  Greek  indices.  We  use  the  notation 


and  assume  that  the  particles  of  each  constituent  all  occupy  the  same 
position  at  time  t  .  We  refer  to  this  position  at  time  t  as  x^ 
and  write 


(2.3) 


The 


velocity  vectors  at  the  point  x^  at  time  t  are  given  by 


«  a 
a  D  xi 
Vi  =  Dt 


(«-  1>2 . v) 


(2.4) 


a 


where  D/Dt  denotes  differentiation  with  respect  to  t  holding 
fixed  in  the  cr  continuum.  This  operator  may  be  written  in  the  form 


i.  =  a  +  v*  -a 

Dt  at  m  ax 


m 


(2-5) 


Acceleration  vectors  at  time  t  are 


®  or  .a  .  a 

j>  _  iTi . 

l  Dt  *  at  +  Vm  ax 

m 


(a-  l|2,...,v) 


(2.6) 


the  mass  densities  at  time  t  are  p  and 

Kor 


m  =  +  -2_  (p  v“)  . 

a  9t  axk  k' 


(2.7) 


We  define  total  density  p  and  mean  velocity  v^  by  the  equations 


p  =  S  p  , 

a=l 


(2.8) 


v  a 

pvi  =  Z  Vi 
1  a-l  a  1 


(2.9) 


We  also  put 


(') . 2U„ ii_it v  au  , 

'  '  Dt  at  m  ax  ’ 


m 


(2.10) 


3. 


and  observe  that 


Of 

_D_ 

Dt 


D_ 

Dt 


v 

+  E 
af=l 


Deformation  gradients  for  each  constituent  are  defined  by 


(2- 12) 


The  residual  energy  equation  for  the  mixture  given  by  Mills  [3)  and 

recast  in  a  slightly  different  form  by  Green  and  Naghdi  [l]  can  be  written 
+ 

as 


pr  -  a 


k,k 


DU  * 
'  P  Dt  •  $ 


-  E  (tt* 
ar=l 


Of 

v.  + 
1 


Of  Of  V 

°ki  Vi,k^  " 


(2.13) 


where  a  comma  denotes  partial  differentiation  with  respect  to  holding 

t  fixed,  r  is  the  heat  supply  function,  q^  the  heat  flux  vector  and 
aki  are  par^ial  stresses.  Also 


U 


v 
■=  E 

Of=l 


u 


v 

»  E 

Of*l 


[p 

3xk  « 


\U0 


\=vk-\ 


(2.3-4) 


where  is  the  internal  energy  per  unit  mass  of  the  constituent  a 
allowing  for  all  interactions  between  this  and  the  other  constituents. 
In  addition 


+Our  equations  of  motion  and  energy  and  our  entropy  inequality  can  be 
shown  to  be  equivalent  to  those  given  by  Truesdell  and  Toupin  [4]  and 
Truesdell  [5!. 


4. 
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(2.15) 


j.  a 

^ki  .or.  Of  l  ,  of  Vv 

^ - +  p  (F.-f.)  =  tt.  +  r  m  (v.-v. ) 

3x,  K<y  i  i'  x  2a'i  i' 


(«  i  v) 


5gki  ,_v  „v.  v  1  v  of  v 

—  .  P^F.-f,)  *  ni  ♦  2  J1"a(''i-Ti) 


where  are  body  forces  per  unit  mass  and 


E  TT%  0  ,  Z  (a“  -  a"  )  -  o  ,  Z  m  -  0 

a-i  1  o-l  K1  o*l  a 


(2.16) 


We  observe  that 


DU  ,  D  Utt  , 

p5t*‘  -J  V  +  ”,  V 


(2.17) 


If  T(>  0)  is  the  temperature  and  entropy  per  unit  mass  of 
the  constituent  of  of  the  mixture,  then  the  entropy  inequality  is  given 
by 

pT  jj§  +  T  Y  -  pr  +  T(^)>k  5  0,  (2.18) 

where 

ps  =  z  p  S  ,  y  =  z  ±r  (p.  \  8)  .  (2.19) 

0=1  a  a  «-l  3xk  ®  *  * 

Let  A^,  A  be  the  free  energy  functions  defined  by 

v 

A  =  U  -  TS  ,  pA  *  Z  p  A  ■  p(U-TS)  .  (2.20) 

o  o  o  O  o 


5. 


Then,  from  (2.13)  and  (2.l8),  we  have 


/DA  _  DT» 

"Pfe  +  S  — ) 


~  *k  9T  "  ,  a  a  a  '  n  .  . 

‘  6  ■TST,.I,("i  V1 +  °ki  sr>  - 0 


k  a=l 


(2.21) 


where 


p„\s„ 


(2.22) 


Also 


DA  V  D  A 

p^r+  ®3  E  (p  -7T7^  +  m  A  ) 
H  Dt  i VKa  Dt  a  a 


(2  -23) 


An  examination  of  (2.21)  suggests  that  we  put 


a  90  ~~ot  ot  ,  ,oi  ~ol 

nk~f^+nk  »  aki  “  8ki  0  +  CTki  5 


(2.24) 


where 


^  (A  .  A  )=  p  (A  -  A)  ,  E  0°=  0  . 

6=1  p  a  p  a  a  a*l 


(2.25) 


In  view  of  (2 .24)-(2 .25) ,  (2. 21)  reduces  to 


/DA  Q  DT\  qk  3T  %  a  -a  9viv  n 

'P(Dt  +  S  5t>  "  T  3^  +a^(TTi  Vi  +  aki  ~  0 


(2.26) 


Moreover,  using  (2.20),  (2.22),  (2.24)  and  (2.25),  the  energy  equation 
(2.13)  becomes 


6. 


/DA 

-  P(«T  + 


-  Of 

„  DT  _  DS»  'J  ,-ot  a  - a  °vi\ 

S  Dt  +  T  S.  i  +  cki  ax.  ^ 

or*l  k 


(2.27) 


We  observe  that  the  parts  of  the  partial  stresses  and  diffusive  forces 

which  depend  on  <t>a  do  not  contribute  to  the  equations  of  motion  (2.15), 

v 

the  total  stress  £  cy“  and  the  energy  equation  (2 .27)  . 

1K1 

Since 


we  see  that  occurs  in  the  basic  equations  of  the  theory  only  in  the 
combinations 


S  ,  S  -  S  (ere  i,2,...,v-l)  ,  (2.28) 

Of  v 

and  of  these  only  S  appears  in  the  entropy  inequality.  Moreover, 

S^-S^  occurs  only  in  the  expression  (2.22)  for  q^  and  does  not 
contribute  to  the  equations  of  motion  and  energy,  or  to  the  partial 
stresses  and  diffusive  forces .  They  must  be  specified  by  constitutive 
equations. 
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3 .  Mixture  of  elastic  continua 


We  restrict  attention  to  the  case  when  the  mass  of  each  constituent 

is  conserved  so  that 


m  =  0 

a 


(3.1) 


Two  methods  seem  to  be  available  for  discussing  constitutive  equations. 

Of  a 

In  the  first  we  assume  that  A  ,  a.  (a*l,..,v),  tt.  (<v-  1, . .  ,v-l) , 

Qf  K1  1 

/  V\  ^ 

(and  hence  tt^),  S  and  are  functions  of 


B  01  *  a 

>  Fij  ,  — g*  (P  «  1, . .  ,v)  , 

3Xk 


(3.2) 


and  linear  functions  of  degree  one  in 


T,k’VVk  • 


(3.3) 


We  might  regard  a  mixture  of  v  elastic  continua  as  one  whose  constitutive 
equations  involve  nonlinear  functions  of  the  quantities  in  (3.3)  as  well 
as  those  in  (3.2),  but  we  restrict  our  attention  here  to  functions  linear 
in  the  variables  (3*3) »  If  we  use  these  constitutive  assumptions  in  the 
inequality  (2 .21)  we  can  show  that  A  and  S  reduce  to  functions  of 


t  ,  f: 


(3.4) 


and  that  restrictions  are  placed  on  A#  and  other  quantities.  In  the 


© 


general  ease  of  v  constituents  there  is  some  algebraic  complexity  in  making 

explicit  deductions  about  the  form  of  from  these  restrictions,  although 

some  information  can  be  obtained  about  the  "equilibrium"  values  of  A  ,  i.e., 

P  v 

when  -  v^  vanish. 

In  the  second  method  we  assume  that  A^  (and  hence  A  and  0a),  S, 
(cr=l,..,v),  q*,  tt?  (or=  1, . .  ,v-l) ,  and  hence  tt^  ,  are  functions  of 
the  variables  in  (3*2)  and  linear  functions  of  degree  one  in  the  variables 

(3.3) -  We  use  the  inequality  (2.26)  to  place  restrictions  on  these  assump¬ 
tions  and  again  we  find  that  A  and  S  reduce  to  functions  of  the  quantities 

(3.4) .  Inspection  of  (2.24)  shows  that  while  will  then  depend  on  the 

functions  (3-2)  and  also  depend  linearly  of  degree  one  on  (3.3),  the 
a 

quantities  are  in  addition  dependent  on 


32T 


3xr  9Xk 


-2J3 
5  FiJ 


0  v 
V  -  v 
i,k  i,k 


(3.5) 


if  we  ignore,  for  the  moment,  any  restrictions  which  might  arise  from 
(2.21).  We  add  the  additional  restriction  that  tt®  is  a  function 
only  of  the  variables  (3.2)  and  a  linear  function  of  degree  one  of  the 
variables  (3-3).  With  this  extra  condition,  coupled  with  the  restrictions 
already  found  on  A  ,  it  can  be  shown  that  A^  (and  hence  Qa)  reduce  to 
functions  of  the  variables  (3.4).  Then  rr®  depend  on  the  same 

variables  (3.2)  and  (3*3)  as  in  the  first  method.  The  remaining  restric¬ 
tion  can  be  found  from  (2.21)  (with  (2.23)),  or  from  (2.26).  We  use  the 
second  method  here  which  is  slightly  more  restrictive  than  the  first,  and 
we  quote  the  final  resultst  Thus,  in  addition  to  the  result  that  A^ 
is  a  function  of  the  variables  (3.4),  we  have 


+ There  is,  however,  no  essential  difficulty  in  adopting  the  first 
procedure . 


9. 


J 


(3.6) 


C_  ;  Po:^_  dA 

"  "a=l  V  3T  '  '  9T 


T®  =  s  p  11.  f* 

ki  P-1  P*  BF*  k<5 
Ij 


(3.7) 


v  M  3F?.  dXP  3A-  9F?.  ax“ 

a  j;  /  _ a  ij  m  ij  tn 

"k  Vl  P“  SF?,  3X?  8?  '  Pfi  8J?,  5X“  3x“ 


4  A 

ij  m 


ij  m 


_  ap  /  P  V\  ,  a  3T 

+  £  c, H  ( v  -  v  )  +  b,  - —  , 

.  .  km  v  m  m  km  Sx 
g*l  m 


(3.8) 


*.D  2Lt£  DJ(vP.v'') 

k  km  9x  a  .  km  m  m' 
m  P=1 


(3.9) 


where 


i  c«*=0 


\  bkm  "  0 

ar=l 


(3.10) 


and  the  coefficients  in  (3-8)  and  (3-9)  are  functions  of  the  quantities 

in  (3.2).  These  coefficients  satisfy  inequalities  arising  from  a 

B  v  /  \ 

quadratic  inequality  in  T  .  and  v''  -  v,  (p  =  1, . . .  ,v-l)  but  these  are 

$  K  K  K 

not  recorded  here.  Finally,  we  note  that  the  above  constitutive  equations 

are  subject  to  the  usual  invariance  conditions  under  superposed  rigid  body 

motions  of  the  whole  mixture.  For  example  A  reduces  to  the  new  form 

Of 


dx“ 

Aa  “  H §  »  T)  (P-l,...,v) 

a  a  ax“  axj 


(3.11) 


10. 


We  observe  that  if  we  had  adopted  the  first  procedure  outlined  at 

the  beginning  of  this  section  the  final  constitutive  equations  would  have 

been  more  complicated.  It  is  still  possible  to  add  further  terms  to  the 

constitutive  equations  of  the  type  discussed  by  Green  and  Naghdi  [6], 

but  these  yield  partial  stress  contributions  which  are  not  necessarily 

derivable  from  the  energy  functions  A  and  are  omitted  here.  Such 

a 

terms  would  influence  the  values  of  the  partial  stresses  and  diffusive 
forces,  but  would  make  no  contribution  to  basic  equations  of  motion  and 
energy  or  to  the  total  stress. 
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